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1. Introduction 

The production of multiple electroweak bosons provides important channels to test LHC data 
against the well-established Standard Model (SM). As the mechanism of electroweak symmetry 
breaking is currently unknown, precise predictions of electroweak gauge boson production rates 
are very important for quantifying deviations from the SM. While the identification of electroweak 
bosons, especially in their leptonic decay channels, is well established experimentally, the theoret- 
ical uncertainties we face at hadron colliders are manifold. 

Total rates and differential distributions suffer from severe theoretical ambiguities if they are 
limited to the leading order approximation in the perturbative series' expansion. The uncertain- 
ties intrinsic to fixed order calculations are conventionally assessed by investigating variations of 
renormalization and factorization scales. The size of these theoretical uncertainties turns out to 
be particularly large when we deal with a LO color singlet final state configuration. This is real- 
ized in pure electroweak production processes such as pp — W or pp — VVV (V = W^,Z, 7). 
In these channels the truncation of the perturbative series at LO means that we force QCD to the 
odd situation where no color gets transferred from the initial to the final state. In other words, the 
next-to-leading order QCD corrections to pp VV,VVV [1-3] are large due to kinematically and 
dynamically unsuppressed initial state parton radiation at large center-of-mass energies. In fact, ex- 
tra parton emission {i.e. the real emission contribution) is required for the cross sections' infrared 
safety only in the soft and collinear limit. Hence, a reasonable question to ask is whether we should 
include hard and resolvable jet emission at all. Vetoing the additional jet activity, i.e. computing 
the exclusive hadronic cross section to next-to-leading order, has been shown to quickly produce 
unreUable results [4-7] and should therefore be considered as a path which should be avoided phe- 
nomenologically. Even worse from a theoretical point of view, an appropriate choice of the jet 
veto scale can be utilized to balance the scale dependencies of the (IR-regulated) real emission 
contribution with the remaining parts of the NLO cross section. For e.g. NLO QCD W 7 production 
processes {y/s = 14 TeV) this is achieved by an experimentally reasonable veto scale choice of 
pf^° ~ 50 GeV. 

However, many search strategies for new physics in the context of multiboson production 
build upon specific final state kinematical configurations and features of the scattering amplitude 
which are special to the leading order approximation only. An example is the search for anomalous 
trilinear couplings in W^7 production [8] via the identification of the so-called radiation zero [9], 
which is a coherence effect subject to distortion when additional parton radiation is taken into 
account [10]. 

Therefore, two serious limitations are seemingly avoided by imposing a jet veto in anomalous 
couplings searches: (i) the cross section seems to be perturbatively stable, and (ii), the radiation 
zero (and the sensitivity to anomalous couplings in total) improves compared to the inclusive NLO 
computation. However, this comes at the price of an uncertain theoretical prediction which can not 
be trusted at the currently known order of the perturbative series' expansion. 

From this perspective the computation of multiboson-f-jet cross sections is important for two 
reasons. On the one hand, adding NLO precision to the one-jet-inclusive production provides an 
important piece of the full NNLO multiboson production. On the other hand, theoretical uncertain- 
ties of the vetoed cross section can be addressed with sufficient precision. As a side benefit we can 
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Figure 1: Sample Feynman graphs contributing to pp — > Wy+jet production at NLO QCD. Left: Bom and 
virtual contributions (the crosses mark possible photon couplings). Right: real emission contribution. The 
pp WZ+jet Feynman graphs follow by replacing the photon with the effective (7*,Z) decay current. The 
figures are taken from Ref. [4]. 



precisely analyze the potential sensitivity to anomalous couplings which arises from recoils against 
important initial state radiation. 

2. Inclusive searches for anomalous couplings 

We consider the most general C and P invariant extensions of the electroweak gauge sector 
[1 1], which modifies the trilinear vertex 

^^^^ = -ie [WI.W^A' - WlAyW^" - KyWlW.F^"' + ^ (2. 1) 

and the trilinear WWZ vertex 

^wwz = -iecot 0,. [gf (W;,W^Z^ - WjZv W^^) - KzW^WyZ^''' + i^w] Wv^Z^^] . (2.2) 

Unitarity of cross sections at high energy scales requires the anomalous parameters ( Ky, Kz,Xy, Xz,gf) 
to be understood as formfactors with a momentum dependence such that for large momentum 
transfers we recover the SM from Eqs. (2.1) and (2.2). Conventional choices are dipole pro- 
files [12], e.g. Ky = I + Ak'o/(1 + m^Yy/ A^), where mwy denotes the mass of the invariant Wy 
system and A relates to the cut-off scale underlying the effective theory of Eq. (2.1). The anoma- 
lous couplings modify a subset of Feynman graphs which contribute to pp ^ £^^^7+ jet + X and 
pp — > i^^f^^i'^ ~^ + ^ ■ The contributing Feynman graphs at NLO QCD are indicated in 
Fig. 1. 

We compute the NLO hadronic cross section by straightforward application of the Catani- 
Seymour dipole subtraction [13]. The loop contributions are evaluated using the Passarino-Veltman 
scheme up to four-point functions [14] and the Denner-Dittmaier reduction [15] for five point 
integrals and we perform various cross checks to validate our implementation (for details see 
Refs. [4-7]). In Fig. 2 we show the total inclusive pp — )• £^Vf7+jet cross section for anomalous 
production. The chosen cuts are adopted from typical searches for anomalous couplings which 
seek to project out the phase space regions with large momentum transfers in the WIV7 (WWZ) 
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Figure 2: Inclusive NLO QCD pp — > £~^^7+jet+X cross section (left, in fb) and K-factor contours (right) 
at the LHC after applying typical search cuts, as a function of anomalous coupling parameters [4, 12, 16]. 
The cutoff scale in the dipole form factor is chosen as A = 2 TeV. Figures taken from Ref. [4]. 



couplings, thus enhancing the sensitivity to anomalous trilinear gauge couplings. In particular, this 
amounts to hard photons {p\ > 100 GeV) or Zs recoiling against the Ws (in the following we focus 
on the vertex [4], for further details on WWZ see Ref. [5]). 

The QCD corrections turn out to be most sizable for anomalous couplings choices close to 
the SM (A; = 0, K", = 1. This is a consequence of the differential K factor being large around the 
Pj threshold while less significant in the distribution's tail. Therefore, large anomalous couplings 
choices, for which the cross section is significantly enhanced in the pr tails, receive a less sizable 
relative correction compared to SM-like choices. The photon's transverse momentum distribution. 
Fig. 3, is typically employed to establish exclusion bounds on anomalous couplings [16]. Even 
though the (rather) inclusive NLO cross sections' scale dependences of ~ 10% turn out to be 
comparable to the relative difference induced by anomalous couplings {cf. Fig. 2), the distributions 




Figure 3: Photon transverse momentum distribution and uncertainty band for SM and anomalous produc- 
tion at NLO QCD from Ref. [4]. Dynamical scale choices have no sizable impact. 
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reveal a significant dependence on the anomalous parameters, especially on A, in Eqs. (2.1), (2.2). 
We perform a binned log-likelihood hypothesis test (taking into account the shape uncertainty that 
arises from QCD scale uncertainty as a nuisance parameter), cf. Ref. [17], to estimate to what extent 
inclusive searches can help to improve the currently pursued analyses strategies ^ We find that 
parameter choices at the edge of the cun^ently allowed LEP bounds on anomalous WVK/ interactions 
can be excluded at the 3a level for ^ ~ 25 fb^^ running the LHC at 14 TeV. Note that this result 
is partly due to the small uncertainty band which follows from the higher NLO precision of the pT 
distributions. 

All results are obtained with an updated version of the Vbfnlo package [19], which is pub- 
licly available. 

3. VFyy+jet production at NLO: Toward anomalous quartic couplings 

Among the triple vector boson production channels, Wyy has turned out to be of particular 
interest. It is the channel with the largest K factor (~ 4) for the integrated cross section, due to the 
radiation zero, and it is also a very promising channel to measure anomalous quaitic couplings [20]. 
The remaining scale uncertainties at NLO QCD [3] are due to unbalanced gluon-induced real ra- 
diation contributions appearing first at NLO, which are computed at LO, e.g., gq — )• Wyy q. The 
observation of the radiation zero is also obscured, similar to production, by additional real QCD 
radiation, W/y-f-jet, as part of the NLO calculation. In Ref. [3] it was shown that an additional jet 
veto-cut might help in the detection of the radiation zero making visible the dip and also reducing 
the scale uncertainties for the relevant distributions. However, this procedure raises the question 
of the reliability of the predictions due to the aforementioned problem with the exclusive vetoed 
samples. To realistically assess the uncertainties, also concerning anomalous coupling searches, 
and as an important step towards a NNLO QCD calculation of Wyy, we have calculated Wyy+jet 
production at NLO QCD. This is the first calculation falling in the category of Wy-njet production 
and includes the evaluation of the complex hexagon virtual amplitudes, which poses a challenge not 
only at the level of the analytical calculation, but also on the level of required CPU time computing 
the full 2 — )• 4 NLO matrix element. 

For the virtual contributions we use the routines computed in Ref. [21] which employ Feyn- 
Calc [22] and FeynArts [23] in an in-house framework. At the numerical evaluation level, we 
split the virtual contributions into fermionic loops (Virtual-fermionbox) and bosonic contributions 
with one (Virtual-box), two (Virtual-pentagons) and three (Virtual-hexagons) electroweak vector 
bosons attached to the quark line. This procedure allows us to drastically reduce the time spent 
in evaluating the part containing hexagon diagrams as explained in Refs. [7,21]. The numerical 
stabihty of the hexagons' diagram evaluation is discussed in detail in Ref. [21]. 

3.1 Results 

For the cross sections and distributions of Tab. 1 and Figs. 4 and 5 we use the CTIO parton 
distribution set [24] with a,(mz) = 0.1 18 at NLO, and the CTEQ6L1 set [25] with as{mz) = 0.130 

'Note that for very hard photons > 100 GeV there is only minimal pollution of misidentified jets and our analysis 
is essentially background-free [18]. 
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Figure 4: Scale variation of the t vyy+i&i production cross sections at the LHC (£ = e, ji) taken from 
Ref. [7]. The cuts are described in the text and we choose = jj.f = mwyy as central dynamical reference 
scale. The right panel shows the individual contributions to the NLO cross section according to our classifi- 
cation of topologies. We also show results where we have applied a veto on events with two identified jets 
having both a transverse momentum larger than 50 GeV. 



W/^77+jet 1.191 1.754 1.47 Tevatron 
lV+77+jet 4.640 6.634 L43 
W-77+jet 3.803 5.644 1.48 

(-) _ _ (-) 

Table 1: Total LO and NLO cross sections and K factors for P P e v^-yy+jet+X and P P 
e+Ve77+jet+X at the Tevatron and at the LHC. Results taken from Ref. [7]. 



at LO. Further details on the parameter choices can be found in Ref. [7]. Again, we consider 
decays to the first two lepton generations, i.e. the decays W — ?■ eVe,/iV^ and these contributions 
have been summed in Fig. 5. 

We choose inclusive cuts to study the impact of the inclusive QCD corrections in a general 
setting: p{. > 20 GeV, p[ > 20 GeV (10 GeV at the Tevatron), > 20 GeV (10 GeV at the 
Tevatron), |TJf |, ITJ7I < 2.5, and an azimuthal angle-pseudorapidity plane separation R(y = {A^f^ + 
Arij^y/^ > 0.4. For the separation of the charged lepton from observable jets, we choose Rij > 0.4 
and we require Ryy > 0.4 for the diphoton separation. Besides the photon-parton isolation criterion 
according to Ref. [26] with the separation parameter 5o = 0.7, we also require a separation between 
photons and identified jets of Ryj > 0.7. 

We compute total K factors of 1.43 (1.48) for W+77+jet (W"77+jet) production at the LHC, 
values which are quite typical for multiboson+jet production as found in Refs. [4,5,27]. This mod- 
erate K factor (as compared to corrections of ~ 400% for 1^77 production) indicates, as expected, 
that the iy77-f-jet production channel is not affected by radiation zero cancellations, since the bulk 
of the W77+jet cross section is due to gluon induced processes which do not sport a radiation zero. 

The scale dependence of the VK+777 and W^Yjj production cross sections turns out to be 
modest: when comparing plr = pLp = ^mwyy for ^ = 0.5 and = 2, we find differences of 10.8% 
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Figure 5: Differential rcaxpj and niy/yy distribution for inclusive and exclusive W 77+jet production. 
Figures taken from Ref. [7]. 



(12.0%), respectively. 

The phase space dependence of the QCD corrections is non-trivial and sizable (we again 
choose Hr = Hf = mwyy)- Vetoed real-emission distributions are plagued with large uncertain- 
ties — a characteristic trait well-known from W-f-jet phenomenology [6,27]. Additional parton 
emission modifies the transverse momentum and invariant mass spectra among other distributions 
such as minimum separations, etc. The leading jet becomes slightly harder at NLO as can be 
inferred from the differential K factor in the bottom panel of Fig. 5. When comparing precisely 
measured distributions in this channel against LO Monte Carlo predictions, the not-included QCD 
corrections could be misinterpreted for anomalous electroweak trilinear or quartic couplings [4,5,8] 
arising from new interactions beyond the SM. 
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